Calculations of spin induced transport in ferromagnets 
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Based on first-principles density functional calculations, a general approach for determining and 
analyzing the degree of spin polarization (P) in ferromagnets is presented. The approach employs 
the so-called tetrahedron method to evaluate the Fermi surface integrations of P in both ballistic 
and diffusive regimes. The validity of the method is examined by comparing the calculated P values 
for Fe and Ni with the experiment. The method is shown to yield highly accurate results with 
minimal computational effort. Within our approach, it is also possible to systematically analyze the 
contributions of various types of electronic states to the spin induced transport. As a case study, the 
transport properties of the soft-ferromagnet CeMnNi4 are investigated in order to explain the origin 
of the existing difference between the experimental and theoretical values of P in this intermetallic 
compound. 

PACS numbers: 75.25.+Z, 85.75.-d, 73.40.Gk 
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INTRODUCTION 

The current upsurge of interest in metal-based as 
well as semiconductor-based spintronics materials mainly 
originates from their increasing applications in the gi- 
ant magnetoresistance (GMR) as well as the tunneling 
magnetoresistance (TMR) based memory devices [3, 0] ■ 
Common to all these devices is the requirement of con- 
trolled transfer or injection of spin-polarized current from 
a ferromagnet into a normal metal Q or a semiconduc- 
tor For such a spin transport to be effective, it 
is mandatory to maintain a relative imbalance between 
spin-up and spin-down electrons that are available at 
Fermi level, ep. By definition, degree of spin polariza- 
tion, P, is described as the extent to which the transport 
current is spin polarized. Depending on the degree and 
asymmetry of spin-splitting, different values of P, rang- 
ing from ~ 40% for a conventional ferromagnetic metal 
such as Fe to ~ 100% for the so-called half-metallic sys- 
tems [H have been reported. In the latter, all the trans- 
port current is carried by either spin-up or spin-down 
electrons. Heusler alloys [6| , chalcopyrites [7|] are the best 
examples of half-metallic ferromagnets whose magnetic 
behaviour as well as electronic structure have been ex- 
tensively studied both experimentally and theoretically. 

In practice, the spectroscopic techniques, such as spin- 
polarized photoemission point contact Andreev re- 
flection (PCAR) (|| and tunnel junction (TJ) [l(J, are 
widely used to probe P at ef of a ferromagnet. However, 
except for the case of half metallic systems, different ex- 
periments may measure different values of P for the same 
candidate material The reason is that, in a realistic 



experiment, the observed tunneling current is dependent 
on both the mean free path of electrons in the ferromag- 
net, I, and the characteristic size of the contact d. If d is 
smaller than the electron mean free path I, the electrons 
flow through the contact ballistically. In the opposite 



case, when d 3> I, they perform a diffusive motion. Con- 
sequently, the obtained values of P, in each regime, can 
be considerably different. In this connection, it is impor- 
tant to compare the measurements with the appropriate 
calculations, mimicking the experiment in question. 

Within the classical Bloch-Boltzmann theory P 
can be defined as (Jj — ■/;)/( Jf + Jj.) where J^m are 
the spin-dependent current densities, passing through the 
contact. Neglecting the state-dependent transmittance 
of the barrier at the point of contact, Mazin [l3| has 
formulated a more practical definition for P as, 



Pn = 



(Nv F ), + (Nv F ) 



(1) 



where, N and vf are the density of states (DOS) at Ep 
and Fermi velocity of electrons with spin a (f or J.), re- 
spectively. Interestingly, all the quantities in eq.[T]are re- 
lated to the band structure of the ferromagnet. Further- 
more, this definition allows direct comparison between 
different experiments and theory. In the simplest case of 
n = 0, the spin polarization calculated from eq. H](Po) is 
in accordance with that measured by spin-resolved pho- 
toemission techniques. Correspondingly, it can be shown 
that, the transport experiments, such as PCAR and TJ, 
measure a higher order of P n (n = 1 or 2) which includes 
the Fermi velocity (see eq. [T]). In the ballistic regime 
(<i<C I), following Sharvin's approach [idj or Landauer- 
Biitikker formalism 15], one can easily prove that the 
conductance of the contact G„ and, consequently J a are 
proportional to (Nv F ) (T [H, El El]. Thus, the ex- 
perimentally observed P in the ballistic regime is ex- 
pected to be equal with P\. On the other hand, us- 
ing Bloch-Boltzmann equation, one can prove that, G CT 
and, consequently J a in the diffusive regime (d ^> I) are 
proportional to ^Nvp') a T a , where r is the relaxation 
time [H 0|. Hence, assuming the same r a in both 
spin channels, the spin polarization observed by trans- 
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port techniques in diffusive regime, is expected to be in 
accordance with P 2 - 

Computationally, it is rather difficult to calculate 
(Nv F ) for n > 0. The problem is due to the presence 
of gradient term up=Vk£, which makes calculations very 
sensitive to the number of A;— points, used for sampling 
the Brillouin zone. Consequently, in most of the the- 
oretical works, it is preferred to compare experimental 
data with the calculated Pq. However, such a compari- 
son may lead to a misleading conclusion in describing the 
transport properties of the material under consideration. 
In this work, we present a simple and general approach, 



based on the so called tetrahedron method [17|, [18[ , for 
evaluating P n in both ballistic and diffusive limits with 
high accuracy. First, the details of the method will be ex- 
plained and then the accuracy of P n values for Fe and Ni 
will be examined by comparing our results with the ex- 
periment. Finally, as a case study, the transport proper- 
ties of the recently discovered soft ferromagnet CeMnNi4 
are investigated, in order to find out the underlying rea- 
sons for the pronounced difference between the experi- 
mental [l9| and theoretical [2^, 21, [23] values of P in 



this particular intermetallic compound. 

METHODOLOGY 

As a general definition, the expectation term of eq. [TJ 
{Nv n ) a , can be expressed as [l3j |. 



KaYl [ v kxJi s k\a - e F )d 3 k 



(2) 



where, £kAo- is the energy of an electron in the band A 
with spin a and the wave vector k. At this point, to 
avoid any confusion, we neglect the subscripts A and a 
and simplify the integration part of eq[2]as, 



3{e F ) 



<ISf 



A(k) 



(3) 



where, A(k) corresponds to v F = |Vk£|"- 

For evaluation of the integral j(sp) over Fermi surface, 
Sf , we consider the following scheme: First we define 
an equispaced grid in reciprocal space on which the lat- 
tice vectors of submesh are obtained by dividing a set of 
primitive reciprocal vectors by the integers m, n? and 
773. Then, utilizing the Monkhorst-Pack method [23j a 
set of irreducible fc— points from {774,77,2, 713} is chosen 
for which the eigenvalues of energy, £k, are calculated. 
For the fc— points, which are not in the irreducible part 
of reciprocal space, symmetry can be used to determine 
their corresponding £k in the irreducible zone. 

Next, each subcell is divided into six tetrahedra of 
same volume. As illustrated in Fig [TJ we choose one 





FIG. 1: Schematic illustration of (a) a subcell and its division 
into six tetrahedrons of same volume and (b) one tetrahedron 
spanned out by the vectors ko, ki, k2, and k3 with the cor- 
responding eigenvalues £0, £1, £2 and £3, respectively. The 
hatched plane is the plane with e = ef- 



main diagonal of a subcell as common edge of all six 
tetrahedra. In order to minimize interpolation distance, 
the shortest main diagonal is chosen. Additionally, each 
tetrahedron is defined by its four corners, one of which 
is, without any restriction of generality, at the origin, 
ko = 0, and the remaining three are at ki [i = 1 — 3). 
With some rearrangement, the four ki are ordered so 
that their calculated eigenvalues £i obey the inequality 
£0 < £1 < £2 < £3- 

Finally, within each tetrahedron, e = ep is interpo- 
lated by a linear function, 



£ = £0 + b • k, 



(4) 



where, b can be described in terms of £i and k^. In this 
order, the triple r^, contragradient to k^ are defined as, 



ri 



k 2 x k 3 
V 



k 3 xki ki x k 2 

r 2 = — - — ; r 3 = — - — (5) 



V 



V 



where V is the volume of subcell, V = |ki • (k 2 x IC3) | , 
and r, satisfy the constraint, r; • kj — Sij. Using eq. b 
can be expressed as, 



Si r; 



(6) 



Back to integral[3l A(k) can be written as a function of 
b, A(k) — |b|". Accordingly, j(£^) over one tetrahedron 
becomes, 



(7) 



and the integration reduces to the following simple form, 



«t(£f) = £(£f) |bf 



(8) 



where £(£_f) is the cross section of the plane e = Ef with 
the tetrahedron (see Fig [JJ . Inserting eq. [5] into eq. [S] 
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and using the shorthand notation e.y for Si — Sj , one can 
easily prove that it(£f) becomes, 



, \ V (ef - £o) 2 ,, 
«t(£f) = — lb| 



if e < e F < £i, 
«t(£_f) = — 



2 £iq£20£30 



£l0£20£30 

if £\ < ep < £2 j and 



£lo£2l£3l 



/ \ V" (£3 — £f ) 2 111 

1 T (£f) = ~n l b l 

^ £3Q£3l£32 



(9) 



(10) 



(11) 



if £2 < £f < £3- 

Here, it is important to point out that, within the 
improved tetrahedron method, proposed by Blochl et 
al. 24| , the linear interpolation in eq. QJJ can be mod- 
ified as, 



it(sf) 



V \h\ n 
■ — ■ — x 

2 £20£30 

£10 + 2(e F - £1) - 



(£20 + £3i)(£F - £i) 2 



£2l£31 



(12) 



Accordingly, throughout this work eq. [12] has been con- 
sidered instead. 

At the last step, summation over it(£f) from con- 
tributing tetrahedra results in the evaluation of integral 
j(sf)- Here, it is crucial to point out that, the de- 
pendence of zt(£f) on the geometry of tetrahedrons for 
n > 0, implies that the summation should be carried out 
over all tetrahedra in the reciprocal space, otherwise a 
misleading error will occur in the evaluation of the inte- 
grals. 

In the following, we apply this prescription to cal- 
culate P n for a number of ferromagnetic materials, us- 
ing their electronic band structure obtained from the 
first-principles calculations. While our method is quite 
general and applicable to any system, we first present 
some benchmark calculations of Pq, Pi and P2 for bcc- 
Fe (a= 2.866 A) and fcc-Ni (a= 3.524 A), and compare 
the results with the experimental data, available in the 
literature 0, [i"o| . Then, we compare the P n values for 
CeMnNi4 obtained from the current method with those 
reported in previous theoretical investigations [HI, [i^] as 
well as experiment jl9| . Here, we focus on the transport 
spin-polarization of cubic CeMnNi4 (a= 6.987 A) con- 
taining 4 formula units i.e. 4 Ce, 4 Mn and 16 Ni atoms. 
For the generated irreducible k~ points (see discussion 
above) the eigenvalues, £kAo-, of all the available energy 
bands, A, at Bp are calculated in both spin channel a and 
used for evaluating P n . All the first-principles calcula- 
tions are carried out using the local spin density approx- 
imation (LSDA) method with projector-augmented wave 



TABLE I: Comparison of P values (in %) for Fe, Ni, and 
CeMnNi4, as obtained from our calculations (P„ (n = — 2)), 
and as obtained from PCAR and TJ measurements (Pc and 
Pt, respectively). The table also includes the calculated aver- 
aged Fermi velocities {vf) in both spin channels (in 10 7 cm/s). 



Structure 




Present work 




Experiment 


Po 


Pi P2 


(wf) t 


(v F ) l 


Pc a Pr b 


Fe 


54 


42 37 


3.1 


4.4 


43 40 


Ni 


-80 


-47 15 


5.5 


1.6 


±46.5 23 


CeMnNi 4 


-22 


-5 10 


2.2 


1.5 


±66 c 



"from Ref . © 
6 from Ref. [lj 
c from Ref. 12 



(PAW) potential, as implemented in the VASP code [25( . 
Details of the electronic structure calculations have been 
published elsewhere [20j . 



RESULTS AND DISCUSSION 

Table [TJ summarizes the P values for Fe, Ni, and 
CeMnNi4, as obtained from our calculations (P n , n = 
— 2), and as obtained from PCAR and TJ measure- 
ments (Pc and Pp, respectively). Here, we first con- 
centrate on the transport properties of Fe and Ni. For 
both bulk ferromagnetic metals, the table indicates an 
excellent agreement between Pi (P 2 ) and Pc (Pt) val- 
ues. Such an agreement indicates that, the experimental 
measurements of Pc and Pp for both metals have been 
carried out in ballistic and diffusive regimes, respectively. 
Fortunately, the details of the corresponding experiments 
explicitly confirm this statement (9I. Ildj. 

Furthermore, from table [TJ one immediately notices 
two different trend of spin polarizations between Fe and 
Ni. In the former, the differences between Po, Pi and P2 
values are moderate while in the latter, such differences 
are so considerable that leads to a change in the sign of 
P-2- The reason can be explained in terms of the elec- 
tronic band structure of transition metals at the Fermi 
level. In general, the Fermi surface of transition met- 
als can be distinguished in two areas. The first area is 
characterized by the s— electrons which can be regarded 
as highly mobile carriers with non-localized wavefunc- 
tions and high Fermi velocity, vf- In contrast, the sec- 
ond area is made by electrons which are relatively lo- 
calized in partially occupied d— like states with a large 
effective mass and low vf ■ Thus, it is expected that, the 
d— like states dominate N and, consequently, Pq, while 
the "light" s— electrons contribute substantially to the 
averaged values of Fermi velocity, (vp} a . As a result, Pi 
and P2 may change dramatically if such a contribution in 
one spin channel is much larger than that in the opposite 
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FIG. 2: (Color online) Spin polarized (a) total and 
s— projected DOS, N and N 3 , respectively, as well as (b) 
Sharvin conductance (Nv) and plasma frequency (Nv 2 ') in 
Fe. The up- and down-states are shown by arrow. 



FIG. 3: (Color online) Spin polarized (a) total and 
s— projected DOS, N and N s , respectively, as well as (b) 
Sharvin conductance (Nv) and plasma frequency (^Nv 2 ^ in 
Ni. 



one (e.g. if (uf) t > (vf) l )- 

For a better understanding, we illustrate in figures [5] 
and [31 both spin polarized total and s— projected DOS, 
N a and N*, respectively, as well as the Sharvin conduc- 
tance, (Nv) a , and the plasma frequency, (Aw 2 )^ for Fe 
and Ni , respectively. In Fe, it is evident that, both 
N a and N* are relatively larger than their corresponding 
spin-down values. Thus, the obtained spin polarizations 
in all limits are positive [P n > 0). Additionally, figure^a 
indicates a large s — d hybridization in both spin channels, 
implying that the averaged Fermi velocities, (wf)| and 
(vp), should not be very different. Hence, the weighting 
factors, vp and v F , in eq. [5] do not considerably change 
Pi and P2 in comparison with Po, as noted in table [J and 
can be seen in figure [2}b. 

On the other hand, Ni behaves differently. In this 
case, the spin-up d— bands are fully occupied and deeply 
located below sp, so that there is almost no spin-up 
d— state at ep, whereas, the spin-down d— bands are par- 
tially occupied and have a significant contribution to the 
Ni (see figure O-a) ■ This leads to a large negative po- 
larization of DOS, (P < 0). Here again TVf > N[. 
However, unlike in Fe, the difference between the aver- 
aged Fermi velocities is so considerable that according 
to table [H (u_f}| is almost four times larger than (vp)^. 
This can be attributed to the fact that in the absence 
of "heavy" d— electrons, Nf is mainly s— like and, hence, 
the transport carriers are expected to be very mobile with 
high vf, while in the spin-down channel, the carriers are 
dominantly d— like and almost immobile in transport. In 
the other word, although the number of electrons at ep in 
spin-up channel is much smaller than that in spin-down 



channel (iVf <C iVjJ but they have more tendency to con- 
tribute to the transport [ivp)^ 3> (vp),). Consequently, 
such an effect results in a significant change in Pi and P2 
values (see figure [3]-b). Due to the same reason, the sign 
reversal in P 2 is also expected. 

Having confirmed the accuracy of our method, we next 
investigate the transport properties of CeMnNi4. Both 
experiment and theory [2(J, 21, 22 1 indicate that, this 
intermetallic compound is a soft ferromagnet with a large 
magnetic moment (~ 4.9/xg/Mn.) and reasonably high 
curie temperature (~ 150K). However, there is a signifi- 
cant difference in the experimental and theoretical values 
of P. PCAR measurements by Singh et al. [l9( yield a rel- 
atively large spin polarization, Pc = 66%, whereas the 
previous first-principles calculations, using both LSDA 
and LSDA+U methods, have estimated much lower val- 



ues for P, varying from -21% to 3% [2jJ, |22j. Addition- 
ally, our calculations on this particular material reveal a 
metallic nature with a very low degree of spin polariza- 
tion. According to table H both P x (-5%) and P 2 (10%) 
are extremely far from the experimental data. This im- 
plies that, it is not possible to realize, through the calcu- 
lations, whether the Pc measurements have been carried 
out in ballistic or diffusive regimes. However, It may 
also imply that such a high degree of spin polarization is 
unlikely to be achieved in pure CeMnNi4. 

From the electronic structure point of view, the small- 
ness of P n values is due to an insignificant difference be- 
tween N-i and N± as well as (vp)+ and (vp),. To elucidate 
this situation, we have shown both total and partial DOS 
of CeMnNi4 around ep in Figure 0J The figure clearly 
indicates that Aj is slightly smaller than N± and, hence, 
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FIG. 4: (Color online) Comparison of (a) spin polarized total, 
s-, d- and /-projected DOS, N, N s (scaled by 10), N d 
and , respectively, and (b) calculated Po, Pi and P2 in 
CeMnNi 4 . 

the calculated Po turns out to be a small negative value. 
Also, it is evident from the corresponding partial DOS, 
N d and , that the d— and /—like states dominate N in 
both spin channels. Consequently, most of the electrons 
available at £p, either d— or /—ones, are "Zieawj/" and can 
not effectively contribute to the spin-induced transport. 
Thus, the averaged Fermi velocities (vp) a are expected 
to be very small and close to each other, as denoted in 
table HI Interestingly, our calculations reveal that both 
(vf)« and (vf) ^ are even relatively (considerably) lower 
than the corresponding values in Fe (Ni). The effects, 
described above, lead to this conclusion that the pure 
CeMnNi4 is more likely a semi-metal with poor spin in- 
duced transport properties. That is why the calculated 
Po, Pi and P2 are so small and somehow negligible. Here, 
it is worth mentioning that, the slight difference between 
(uf). and (vf) , in CeMnNi4 are due to the fact that the 
contribution of s— like states to is slightly larger than 
that to N± (see figure |3J-a). 

It is important to point out that, although the ob- 
tained P n values at ep are substantially smaller than the 
experimental data, but for changes of the Fermi level 
by approximately 0.1 eV, values close to experimentally 
ones for Po and Pi are reached (see figure HJ-b). Ex- 
perimentally, this situation can occur if there is a small 
stoichiometric variation in the sample. For example, as 
noted in Ref. [U if during the process of preparation, 
small amount of Mn (e. g. ~ 5%) is replaced by Ni 
one may expect a considerable change in the electronic 
structure and hence the transport properties of the chem- 
ically disordered structure. In the other word, Since the 



exact experimental conditions of synthesis of CeMnNIt 
is still not completely understood, it might be possible 
to attribute the measured Pc value to a sample, con- 
taining small amount of impurities, defects or chemical 
disorders. Thus, we suggest further studies on the role of 
stoichiometric variations on the transport properties of 
CeMnNi 4 . 

CONCLUSION 

A general and accurate scheme based on tetrahedron 
method for determining the degree of spin polarization 
in ferromagnets has been presented. Our approach has 
successfully shown to yield highly accurate spin polariza- 
tion results both in ballistic and diffusive regimes, with 
the least number of fc-points and the modest computa- 
tional efforts. The benchmark calculations on bulk Fe 
and Ni have proven the efficiency of the method in re- 
producing the experimental data in different regimes. It 
turned out that the spin polarization in a ferromagnet 
is dependent on both the number and the type of trans- 
port carriers in each spin channel. That is, the more 
the contribution of s— electrons to the density of states 
in one spin channel, the higher is the degree of spin po- 
larization in both ballistic and diffusive regimes. Finally, 
in agreement with previous first-principles calculations, 
the obtained spin polarization in CeMnNi4 indicated a 
remarkable difference with the experimental data. The 
reason was attributed to the strong d—f hybridization in 
both spin channels, leading to a reduction in contribution 
of s— electrons to the spin induced current. It was sug- 
gested that the possible presence of off-stoichiometery in 
the sample might be responsible for such a discrepancy. 
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